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Abstract 

We discuss the hydrodynamic boundary condition for a superfluid moving tangentially to a rough surface. 
Specifically, we argue that the scattering of quantum fluctuations off surface roughness affects the nature 
of the boundary condition, and that this has important consequences including a theorized critical speed 
and the presence of normal fluid at any nonzero speed, even if the boundary is held at zero temperature. 
This hydrodynamic boundary condition is relevant not only for superfluid helium experiments but also for 
experiments with trapped dilute Bose-Einstein condensates, in particular those involving atomic waveguides 
near surfaces. 

PACS numbers: 
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I. INTRODUCTION 



Much of what is accepted about the nature of superfluid flow and its instability is built on a few 
concepts laid down by Landau'. However, many aspects of superfluidity remain mysterious. One 
such is the discrepancy between Landau's predicted critical velocity in superfluid Helium and that 
observed in experiment^. 

Allen and Misener^, in their letter reporting one of the first observations of superfluidity in 
liquid Helium, suggested that superfluidity is due to a lack of friction of the fluid with the walls^. 
Given that the connection between the nature of superfluidity and its boundary conditions was 
clearly made so early on, it is somewhat surprising how little attention has been paid to super- 
fluid boundary conditions or, more specifically, the boundary conditions of a superflow moving 
tangentially to a rough wall (see however-). 

We believe that these boundary conditions merit closer study. This paper revisits the broadly 
assumed boundary conditions on superflow, and hopefully sheds some light on enduring enigmas 
like the value of the critical velocity in superfluid Helium experiments. In particular, this paper 
looks into how the presence of quantum fluctuations affect the boundary conditions on a superfluid, 
and proposes a resulting effective boundary condition (eq. ©). The motivation for our view that 
quantum fluctuations may significantly alter the boundary conditions stems from our previous cal- 
culation of a drag force on a localized object; we showed the critical velocity effectively becomes 
zero for an infinitely extended geometry due to the quantum fluctuations in a superfluid^^^. This 
discussion should be important for recent experiments directing Bose-condensed atoms through 
waveguides near a surface (see e.g.-^^^i^). It might also provide an explanation for the sensitivity 
to surface roughness seen in a recent experiment where a vibrating grid is immersed in superfluid 
Helium^. Other experimental scenarios involving dilute Bose-Einstein condensates (BECs) and 
superfluid Helium will be discussed further at the end of this article. 

In order for the findings to have broad relevance for current superfluid Helium and dilute BEC 
experiments, the geometry considered in this paper will be that of the superfluid system contained 
by rough walls moving at a certain speed Vg relative to the superfluid component. Wherever it 
is potentially ambiguous, we will refer to the coherent component as the superfluid component, 
the remaining fluid (i.e. excitations) as the normal fluid component, and both components of the 
fluid together as the superfluid, or the superfluid system. The container is assumed to be a finite 
system, such as a spinning torus, that locally resembles a long cylinder. However, the condition of 
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uniformity of the phase in such a geometry does not play a crucial role in the ensuing discussion. 

n. LANDAU'S TWO-FLUID PICTURE 

Let us begin with Landau's conception of superfluidity. Landau^ wrote 'Physically it is obvious 
that the interaction of the fluid with the moving walls of a tube cannot cause a motion of the whole 
liquid at once. The motion must begin with the excitation of the inner movements in the layers of 
liquid close to the walls of the tube, i.e. with the excitation of rotons and phonons in the liquid'. In 
his phenomenological picture, excitation of rotons and phonons was the (sole) mechanism through 
which the walls would interact with the superfluid. Because such excitations could only be created 
if the superfluid's motion relative to the wall was above a certain velocity, this led him to postulate 
the existence of a critical velocity below which the superfluid would flow completely without 
energy dissipation. His reflections on superfluid systems implied the following two conditions at 
the boundary: 

Vn = (1) 

n-Vs = 0, (2) 

where the first is the usual non-slip boundary condition for a viscous fluid, which determines 
the velocity of the normal fluid component composed of the excitations of the fluid; and the 
second is the boundary condition of a perfect potential flow, obeyed by the superfluid (or coherent) 
component with velocity Vg. Khalatnikov later generalized these boundary conditions by including 
the thermal conductivity of the solid boundary'^'^. In this paper, however, we will assume that the 
boundary is perfectly insulating — with zero thermal conductivity — and thus we will not consider 
heat flux across the boundary. 

Landau's critical-velocity argument involves representing the moving 'ground state' of a su- 
perfluid flowing tangentially to the container wall as a straightforward Galilean transform of the 
fluid's true ground state in the non-moving reference frame. Landau makes the implicit assump- 
tion that one can ignore the containing walls' roughness when performing this Galilean transform. 
For if one took into account the boundary roughness then the ground state boosted to give the fluid 
a velocity tangential to the wall is not a possible eigenstate of the system. 

To illustrate, consider a superflow moving in the x-direction. In the absence of a containing 
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wall, the fluid's wavefunction would have the property of Galilean invariance. That is, if 

*o(rl,r2,rl,...)e-*— (3) 



is a solution of the Schrodinger equation with a well defined eigenvalue Eq for the energy, where 
7^,. . . are the position of all the particles, then 

^0(^1 - vt, r2 - vt, - vt, . . .)e ^ n e ^ n (4) 

is also a solution, m being the mass of the particles. This remains true in the presence of a flat 
wall, which can be represented by a potential U (z). If the wall were rough however, it would be 
represented by a potential U {x, y, z) where x and y are coordinate axes parallel to the mean plane 
of the wall. Any eigenfunction of the motionless ground state of the system would be a function of 
all the coordinates of all the particles. Therefore the boosted Galilean state cannot be an eigenstate 
of the system for any non-zero value of the velocity. 

Further, Landau's argumentation does not take into account quantum fluctuations. We will 
show that by including the effects of surface roughness - a feature of all physical walls - and the 
effect of quantum fluctuations, the boundary conditions become much more complicated. 

Let us consider that the walls containing the superfluid system are rough. Then, as mentioned 
above, the state of the flowing superfluid cannot be a boosted ground state: there must be some 
changes in what Landau calls the inner motions close to the walls at any value of the speed, if only 
because disallowing a Galilean boost means the loss of a simple solution. However, although the 
state resulting from a Galilean boost cannot be an eigenstate, this does not necessarily mean that 
there is no eigenstate. As we shall see, the physical picture is quite complex. 

One might attempt to adapt the smooth- wall argument to the rough-wall situation by including 
some boundary layer close to the wall, which would merge continuously with the boosted ground 
state far from the wall. This cannot be the case for several reasons. The most obvious (but not 
direct) argument against this is the following: Think of a classical gas flowing tangentially at a low 
speed (low Reynolds number) with respect to a plane wall; if this wall were perfectly smooth down 
to molecular scales and if the fluid particles were elastically scattered by this wall, there would 
be no exchange of longitudinal momentum between the wall and the classical flowing gas; but if 
the wall were rough, there would necessarily be friction because the perturbation brought to the 
velocity distribution in the gas is not localized in the neighborhood of the wall. The wall friction, 
at least in a permanent state, generates a constant shear stress on the fluid at any arbitrary distance 
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from the wall of the finite system such as the long cylinder considered here (we are ignoring 
transient effects and assume the system has reached its steady state). In this example, no boundary 
layer takes care of the roughness to allow merging with a state of uniform speed far from the wall. 

Returning to the superfluid case, a way of understanding what happens near a rough wall is to 
look at the dilute Bose gas limit, which is valid for recent experiments in ultracold trapped atomic 
gases^^. While this limit is not strictly valid for superfluid liquid Helium, there has been a long 
history of its use to gain insight into superfluid Helium behavior. In this limit an expansion in 
terms of the small parameter in the theory {^/noa^, where no is the number density and a is the 
positive 2-body scattering length) allows one to apply perturbation theory. In this limit, the leading 
order is given by the time-independent Gross-Pitaevskii equation (GPE), a classical equation for 
the amplitude of the mean-field wavefunction \E'(°)(r): 

(f + f/(f) - ^)^(°)(r) + ^|^(°)(f)|2^(°)(^) = 0, (5) 

where /i is the chemical potential, T = —^^^ + ^^^s^ + \'mvs'^, g = ATrh?a/2m, and m is the 
mass of the atoms. A stationary potential in a moving flow would thus experience a force 

Fgpe = -J d^f\¥'\f^\'VU{f). (6) 

Eq. ([5]) reduces in the hydrodynamic limit to the familiar Bernoulli equations of a compressible 
fluid (if one ignores the quantum-pressure term) if the maximum flow velocity does not exceed the 
speed of sound — the critical speed identified by Landau. There is a well defined steady solution 
of the flow equation typically with a thin boundary layer compensating for the rough wall. If the 
maximum flow velocity exceeds Mach 1, the steady solution of the GPE disappears and vortices 
are released^2i2i. The behavior of the GPE agrees almost perfectly with Landau's conception of 
superfluid flow, where no irreversible losses occur below this threshold speed. However the GPE 
describes only part of the picture; it leaves out the zero-point quantum fluctuations. 

III. EFFECT OF QUANTUM FLUCTUATIONS 

In the dilute limit, quantum fluctuations appear at the next order in the expansion parameter 
(y/nooP). (See Appendix for a review of the effect of scattering quantum fluctuations in the dilute 
limit and an outline of the calculation of the resulting drag force on a localized stationary object 
immersed in a superfluid moving at arbitrarily low velocities.) 
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Let us leave rough walls aside for a moment and focus on the effect of quantum fluctuations. 
In a trivial sense, the fluid motion does not induce any transition in the occupation number of the 
quantum state. But one has to define exactly what is meant by occupation number here and, more 
specifically, what state one is referring to. The ground state refers to the state that is immobile 
in the lab frame. In this state, the macroscopic wavefunction corresponds to a uniform state, and 
zero-point fluctuations (defined as the ground state of occupation for certain modes of fluctuations) 
also have a well defined amplitude. Practically, each mode of fluctuation is associated with a 
certain wavenumber. Because of the global Galilean invariance of the system (ground state and 
fluctuations) any fluid motion can also be seen as changing the quantum state of the system by a 
Galilean boost. This Galilean boost defines exactly the quantum state. 

If one assumes the wall to be flat down to molecular scales then, regardless of whether the fluid 
is a classical one at finite temperature or a superfluid at zero temperature, there would be no rub- 
bing between the flowing fluid and the wall itself because no force is generated in a given direction 
if the potential does not vary in that direction. Therefore any tangential stress generated from the 
scattering of quantum fluctuations by the boundary would have to come from a modulation of the 
potential wall normal to the flow direction. 

We now reintroduce rough walls. The eigenstates of this system where quantum fluctuations 
scatter off the rough walls differ from those of a homogeneous Galilean-boosted system of ground 
state plus excitations. It is nevertheless instructive to map the true eigenstates of our system with 
rough walls onto those of the latter system since at large hydrodynamic scales the only non-entropy 
producing states are those of a Galilean-boosted homogeneous superfluid component plus excita- 
tions. The ground eigenstate of our scattering system maps onto the Galilean-boosted ground 
state plus similarly Galilean-boosted excitations, i.e. the normal fluid component, because, in 
general, mapping a ground state of one system onto another system whose Hamiltonian has dif- 
ferent boundary conditions will result in nonzero amplitudes in the excited states of eigenstates 
describing the target system. 

This mechanism for creating excitations differs from Landau's theory of macroscopic excitation 
of quasiparticles that does not depend on the scattering conditions of the rough walls and that re- 
quires speeds larger than Landau's critical speed. The interaction between the rough wall potential 
and the quantum fluctuations is a scattering problem as described in'^ '^, and the relative differ- 
ence of speed between the superfluid and the wall changes the scattering condition. (To discuss 
equilibrium, one must at least implicitly assume that the higher order terms where the excitations 
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interact, i.e. the Landau and Beliaev damping rates (se o ^ • ^ ' ^ and references therein), are 
the source of collisions required for the system to equilibrate.) The temperature difference between 
this normal fluid component and the boundary must depend on the tangential superfluid velocity 
and vanish if this velocity is zero. Below we present some ideas for representing this boundary 
condition. 

IV. EFFECTIVE BOUNDARY CONDITION 
A. Leading order (assuming Vs is constant) 

Taking into account the scattering of quantum fluctuations as described above, the simplest rep- 
resentation of the boundary condition at the interface between a superfluid and its containing wall 
begins with the assumption that the scattering of the quantum fluctuations (also boosted by this lo- 
cal speed) generates a temperature difference between the wall and the superfluid system, resulting 
in some normal fluid component at rest along the boundary even if the boundary is held at T5 = 0. 
This temperature difference is proportional to the square of Vg — the tangential superflow speed 
— since Vg and (— fs) have the same effect. Ultimately, in equilibrium there should be no ongoing 
generation of entropy (we expand upon this idea below). To make a boundary condition out of 
this assumption one has to assume that, in the finite system, if the temperature difference is not 
that at which this equilibrium condition is reached, the tangential speed should induce some net 
transformation of superfluid component into normal fluid component or vice versa at the bound- 
ary. (This is reminiscent of what occurs when there is a normal heat flux at the boundary, often 
referred to as Kapitza resistance^ although it is important to stress that we only consider a per- 
fectly insulating boundary condition in this paper.) One can then write a similar phenomenological 
non-equilibrium fluid mechanics boundary condition (valid for scales much larger than that of the 
surface roughness) for such a situation, yielding the main result of this paper: 

n-{3n-]s) = Oi{T,-T)+(3vl (7) 

where n is the unit vector normal to the average (macroscopic) surface, j„ and js are the mass 
fluxes for the normal fluid component and superfluid component respectively, is the tempera- 
ture at the boundary, and T is the local temperature of the fluid, a and (3 are phenomenological 
parameters which depend on the surface roughness. We arrived at this boundary condition by 
speculating that the equilibrium temperature is proportional to the superfluid component speed vl 
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multiplied by a parameter characterizing the surface roughness and by assuming that Vs is approx- 
imately constant. This tangential boundary condition is to be supplemented by the total normal 
mass flux condition for a solid wall 

n-{]n+3s)=Q (8) 

as well as by the normal condition for energy transport across the boundary'^, which is responsible 
for second sound but has no bearing on the present discussion as we do not consider energy flux 
across the boundary. 

At thermal equilibrium, n ■ jn = (actually j„ = 0) because there should be no continuous 
generation or absorption of entropy. Therefore, n-js = (from eq. ([8])) together with eq. ([7]) yield 
the following expression for equilibrium temperature T of the fluid: 

T = Tk + ^. (9) 

a 

Let us note here the peculiar thermal equilibrium in the situation considered in this paper. If 
Vs = in the situation considered in this paper, thermal equilibrium would imply that the normal 
fluid component (and thus the superfluid system) would be at the same temperature as the bound- 
ary; but at a nonzero Vg the scattering of quantum fluctuations will necessarily force a temperature 
difference between the boundary and superfluid system at equilibrium. That is to say, T would be 
greater than Th in thermal equilibrium. This can occur because the usual argument for uniform T 
at equilibrium as a result of energy exchange in thermodynamics does not apply here, since energy 
is also stored in the motion of the superfluid component. The scattering of quantum fluctuations in 
a moving superfluid would cause conversion from superfluid component into normal fluid compo- 
nent and, with that, the kinetic energy of the superfluid component is converted into heat. By the 
relation eq. ©, this means that £ must be greater than zero. Thus even if Tf, is held at zero and Vg 
is constrained to be nonzero, T will be nonzero. 

B. Higher order correction (local boundary condition) 

Although to leading order Vg can be treated as constant, the roughness of the surface would 
cause some slight variation in Vg across a given flow structure. Eq. ([7]) does not have the flexibility 
to permit equilibrium between the non-constant Vg and a constant wall temperature, a constant tem- 
perature T inside the fluid and no net transformation between normal and superfluid components 
at the wall (as represented by the left-hand side of equation dV])). 
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To permit a nonuniform Vs, we will consider a modification of eq. (|7]). Recall first that in 
thermal equilibrium T > T;, if Vg is nonzero. Because the system must be at thermal equilibrium, 
to take into account scattering by the rough walls one must assume that, apart from the coherent 
superfluid component, there is some amount of normal fluid component moving with the container. 
This is a property of the system as a whole and cannot be written simply as a local boundary 
condition as was attempted in eq. ([7]). Nevertheless we do expect a local relation between the 
various physical parameters at a given point on the wall. To reconcile the existence of both a 
global thermal equilibrium state and a local boundary condition, we suggest that instead of a 
smooth function of a{Tb — T) + /3v'^, the RHS of eq. Q) is modulated by a Heaviside function, i.e. 

n-0n- 1) = («(T, - T) + (3v^^) H[-a{T, - T) - Pv]), (10) 

where H{x) = 1 for positive x and zero otherwise. With this modification, the RHS of eq. dV]) 
vanishes as soon as a{Tb — T) + becomes positive. This rather unusual boundary condition 
seems to agree with all the constraints of the problem, but constitutes only a small correction to 
eq. (|7]). Eq. (flOl) implies that, for a given constant wall temperature and a given field of wall speed 
Vs, equilibrium will be reached at a temperature T in the fluid such that 

T = Tb + -max{v^,) (11) 
a 

where max(f g) is the system-wide maximum speed at the wall. Note again that even if Tf, = 0, 
the temperature of the fluid will go to a nonzero T for a finite system and nonzero Vg. 

V. DISCUSSION AND CONSEQUENCES OF EFFECTIVE BOUNDARY CONDITION 

It is important to understand the difference between this picture and the friction in Landau's 
theory. Here, the friction with the wall is not an ongoing process, unlike the continuous shedding 
of vortices beyond the critical Landau speed; instead the 'friction' in this picture (namely the 
scattering of quantum fluctuations) brings the system reversibly to a state of equilibrium where 
there is some amount of normal fluid component at rest with respect to the wall. This fits well with 
the widely observed existence of persistent currents, something that would be clearly incompatible 
with friction by continuous scattering of quantum fluctuations. This analysis only applies to a 
closed geometry where multiple scattering occurs, like the moving container we alluded to before. 
Indeed all experiments are done in finite containers. However, we expect that equilibrium will be 
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hard to reach in geometries approaching, for example, the classical Stokes situation of a moving 
sphere in a fluid at rest at infinity. In such a case, equilibrium will be reached only when the 
temperature due to the fluid motion near the sphere equilibrates with the temperature 'at infinity', 
at the wall of the container. In a semi-infinite geometry - a superfluid flow over a single wall 
- all of the superfluid component will eventually be converted to the normal fluid component as 
equilibrium will never be reached. This can be likened to the classical fluid mechanism of a flow 
in a semi-infinite geometry — the boundary layer will grow forever. 

There are some important consequences of the boundary condition postulated in the present 
article, eq. dV]). (The correction resulting from the introduction of the Heaviside function (eq. [TO)) 
is slight and has no material effect on the conclusions drawn below.) Firstly, if there is a super- 
fluid flow along a rough boundary, there will always be a small amount of normal fluid component 
at rest with respect to the wall, even if Tb is kept at 0. This normal fluid component, at least in 
principle, could be observed in (immersed) torsion pendulum experiments^^. It would appear as a 
nonlinear correction to the pendulum equation; at increased angular velocity more normal fluid is 
carried by the oscillating discs, so the effective moment of inertia depends on the angular velocity. 
Secondly, this boundary condition provides a critical speed in bulk flow experiments, one above 
which the equilibrium temperature in the flow is larger than the thermodynamic temperature of 
transition. This might provide an alternative to the problematic vortex- shedding hypothesis^ used 
to explain the infamous discrepancy (two to four orders of magnitude difference) between the ob- 
served critical speed and that predicted by Landau. Thirdly, this boundary condition tells us that 
when one slows down the moving boundary, excitations are converted into superfluid component 
(much like superfluid component creation when heat is removed from the boundary in counterflow 
experiments). This newly created superfluid (coherent) component is likely to be not truly coher- 
ent, much like a Bose-Einstein condensate growing out of a nonequilibrium gas, which has a phase 
coherence length growing continuously with time after its inception - typically as (^) from the 
time the condensate is created. In the situation considered, t is the approximate time it takes for 
the newly formed superfluid component to travel a distance ^ from the wall. A rough assumption 
of the phase coherence length, which is approximately the average distance between vortices, for 
our system can be had by substituting t ~ ^/ Us""' for t, where Ug is the component of the superfluid 
velocity normal to the boundary, i.e. Us = n ■ vl . (This substitution requires the assumption that 
the long range order of the phase depends on vortices carried by the superfluid.) Therefore the 

1 /2 

coherence length parallel to the wall is ( j . Almost all vortices will have disappeared if this 
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phase coherence length is larger than the size of the system, and vortices will be everywhere when 
this length is much smaller than the length of the system, with crossover length scale for the size of 
the system being Because ^ has the physical dimensions of a shear viscosity, this condition 
can be likened to a condition on a Reynolds number. 

Although the boundary condition was motivated by superfluid Helium experiments, the 
exquisite control over the atomic interactions and confining potential in experiments with current 
dilute Bose-Einstein condensates make this an ideal medium to test systematically the predictions 
presented in this paper and shed light on superfluid Helium flow experiments past a rough wall. In 
addition to previously mentioned atomic chip experiments where surface roughness is important, a 
possible example would be to drag a laser sheet, where the surface roughness can be precisely con- 
trolled, through a condensate and measure the resulting heating of the gas. Similarily, one could 
modify a version of a previous experiment where a laser was dragged through a condensate^: one 
should be able to give the laser beam roughness using an acoustic-optical modulator and twist the 
beam in the condensate so that it acts like a rough corrugated cylinder rotating in a dilute BEC. 
Another possible experiment would be a superfluid flow over a disordered optical lattice^^-^^iM 
or optical speckle potentiaP^- which could be adjusted and mimic the effects of surface rough- 
ness. Finally, experiments to observe persistent currents in dilute BECs are just beginning to be 
realized^^, and one could, in principle, precisely modify the trap in such a way as to control 
the surface roughness, unlike the experiments on superfluid Helium, and detect the normal fluid 
component of the gas by measuring the temperature and mass flux. 

D.C.R. gratefully acknowledge stimulating discussions with Brian Anderson and Kevin Hen- 
derson. 
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APPENDIX A: SCATTERING OF QUANTUM FLUCTUATIONS 



In the dilute limit, quantum fluctuations appear at the next order in the expansion parameter 
(Vnoo^) defined by expanding the exact bosonic field operator in terms of a coherent field plus a 
small quantum fluctuating operator, i.e. 

^(f,t) = ^«(f,t) + 0(f,t), (Al) 

where "^^^^ is the wavefunction of the condensate affected by the action of the quantum fluctua- 
tions. A stationary potential in a moving superfluid in this dilute limit would then experience a 
force 

^ = - J d^r{\¥^\f)\^ + (0^(f)0(f))r=o)Vt/(f). (A2) 

can be expressed as the sum over excited states of quasiparticle operators and al weighted 
by quasiparticle amplitudes Uk{r) and ffc(r) such that 

m = {Mrl&k - vlif^al) . (A3) 

k' 

Uk{r) and Vk{r) are assumed to satisfy the dynamic equations that, in this dilute limit, govern the 
quantum fluctuations: the Bogoliubov-de Gennes equations (for derivations and explanations see 
for example^^iSi^)^ 

Cukif^ - gi^^'^^^Vkif^ = EkU^ir) (A4) 

tvu{^ - gi^'^^ fuuir) = -EkVkif), (A5) 

where is the eigenvalue associated with momentum state k, C = T + U{f) — fi + 2g\'$^^''\'^, and 
* signifies complex conjugate. The weakly interacting Hamiltonian would then be diagonalized by 
ak and al- In this limit, the behavior of the coherent field modified by the quantum fluctuations is 
given by the Generalized GPB^^di 

{f+U{f)-i^)¥^\7^+g\¥^\f^\''¥'\f^+J2'^9\Mr}\''^^'\r^ = 0. 

k' k' 

(A6) 

The term g J2k' Uk{r)vl{r)^^^^ (r) requires renormalization as the contact potential approximation 
causes it to be ultraviolet divergent. In the final term on the LHS, /(r) = Cfcf^(r), where = 
g J (i^r1\I>(^)(r)p(\[^(^)*(f)Mfc(^) + "^^^^rjvkir)). /(r)\&(^)(f) ensures that the excited modes 
and the condensate are orthogonal to each other.-^^ provide explicit solutions of these equations 
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describing an effective scattering problem for the quantum fluctuations around a stationary object 
in a superfluid moving at any arbitrarily small velocity, and demonstrate the existence of a nonzero 
drag force on the object in such a situation. It is important to note that the situation they consider 
is an infinitely extended system, where the kind of thermal equilibrium described in this article 
cannot be reached. 
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